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squares  relationship  between  the  peak  absorption  (A)  and  the  burst  energy  (F.)  is  given  by 
A (dB)  = 0.01 16£  (Sfu*min)0-5555  with  a correlation  coefficient,  R,  of  0.59.  The  peak 

V C 0 t K - rrSr  , * 

DD  | J AN^T  1 1473  EDITION  OF  I NO  V *5  1 5 OBSOL  E T E UNC  LASSIF I F.D 

SFCURlTY  CLASSIFICATION  OF  THIS  PAGE  When  Data  Fnteied) 


H ' 7- 


UNCLASSIFIED 
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absorption  predicted  by  this  formula  is  within  a factor  of  2 of  that  actually  observed  for 
22  of  the  32  principal  PCA  events  (69  percent)  and  within  a factor  of  3 for  28  of  the  32 
(88  percent).  A correction  for  the  effect  of  the  interplanetary  magnetic  field  on  particle 
propagation  is  considered  in  the  appendix. 
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OBJ  ECTIVE 


Develop  earth  environment-disturbance  forecasting  techniques  to  predict  degradation 
effects  on  military  communications,  navigation,  and  surveillance  and  on  weapon  systems 
which  use  electromagnetic  radiation. 


RESULTS 

Forty-seven  Polar  Cap  Absorption  Events  which  occurred  during  the  nineteenth  and 
twentieth  solar  cycles  were  studied.  Parent-Hare  2.8-GHz  event  energy  was  correlated  with 
peak  riometer  absorption  for  these  events.  Thirty-two  of  the  events  had  riometer  absorption 
>2.0  dB.  The  least-squares  relationship  between  the  peak  absorption  (A)  and  the  burst  energy 
(E)  is  given  by 

A(dB)  = 0.01 16  E(Sfu-min)0-5555 

with  a correlation  coefficient  (R)  of  0.59.  The  peak  absorption  predicted  by  this  formula  is 
within  a factor  of  2 of  that  actually  observed  for  22  of  the  32  principal  PC  A events  (69  per- 
cent) and  within  a factor  of  3 for  28  of  32  t.88  percent)  events. 
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INTRODUCTION 


Polar  Cap  Absorption  (PCA)  events  caused  by  low-energy  (*  10  MeV)  protons  ejected 
from  intense  solar  Hares  can  disrupt  high-frequency  communications  in  the  earth’s  polar 
regions  for  periods  of  days.  I lie  spatial  and  temporal  extent  of  the  PCA  event  make  it  of 
particular  concern  to  the  US  Air  Force.  The  time  delay  between  the  optical  and  radio  obser- 
vation of  the  Hare  and  the  arrival  of  Hare  particles  enables  the  solar  forecaster  at  the  Air 
Force  Global  Weather  Central  (AFGWC)  to  notify  customers,  who  use  the  polar-region  iono- 
sphere in  their  communications  networks,  of  the  impending  occurrence  and  expected  inten- 
sity of  a PCA. 

The  spectral  shape  of  the  Type  IV  microwave  radio  emission  associated  with  large 
Hares  has  proven  a useful  tool  for  predicting  whether  or  not  a given  Hare  will  produce  a prin- 
cipal (>  2.0  dB)  riometer  absorption  event  (ref  1 -3).  Several  methods  (ref  4-12)  have  been 
developed  to  forecast  the  intensity  of  a solar  particle  event  based  upon  other  parameters  of 
the  parent-flare  microwave  burst.  These  methods  are  reviewed  in  table  1. 

The  first  four  studies  attempted  correlations  between  burst  parameters  and  time- 
integrated  proton  fluxes.  Straka  and  Barron  (ref  8)  were  the  first  to  consider  correlations 
with  peak  proton  fluxes  and  riometer  absorption.  In  a follow-up  study  based  upon  data  from 
the  range  of  discrete  frequencies  monitored  at  Sagamore  Hill  Observatory.  Straka  concluded 
that  the  integrated  burst  flux  density  (event  energy)  at  the  longer  centimetre  and  decimetre 
wavelengths  was  a useful  predictor  of  the  size  of  PCA  events  (ref  4).  ('room  obtained  a high 
correlation  (0.8)  between  parent-flare  burst  mean  duration  and  peak  10-MeV  proton  flux  for 
a relatively  large  data  sample,  but  this  method  has  yet  to  receive  an  operational  test  (ref  10). 

An  analysis  by  this  author  of  32  solar  particle  events,  detected  by  the  Goddard  Space  Flight 
Center  Experiment  aboard  Explorers  34  and  41.  revealed  a correlation  coefficient  of  only 
* 0.5  between  the  2.8-GHz.  event  mean  duration  and  the  peak  6-|4-MeV  proton  flux.  The 
methods  developed  by  Newell  (ref  1 1 ) and  Castelli  et  al  (ref  1 2)  are  currently  in  use  at  AFGWC. 


1 . Castelli.  JP.  Aarons,  J,  and  Michael.  GA.  “Flux  Density  Measurements  of  Radio  Bursts  of  Proton- 
Producing  Flares  and  Nonproton  Flares.”  Journal  of  Geophysical  Research,  v 72,  p 544| . ll>o7 

2.  AFCRL  Technical  Report  (>8-0104.  Observation  and  Forecasting  of  Solai  Proton  I vents,  by  JP  Castelli. 
1468 

3.  AFCRL  Technical  Report  70-0425.  The  Prediction  of  Solar  Proton  I vents  Based  on  Solai  Radio  I mis 
sums,  by  W!  O'Brien,  1470 

4.  NASA  Technical  Report  SP-50.  A Review  of  Solar  Cosmic  Ray  Events,  by  WR  Webber.  A.AS  NASA 
Symposium  on  the  Physics  of  Solar  Flares,  edited  by  WN  I less,  p 215.  1 463 

5.  Fletcher.  J 1 > . "Solar  Radio  Emission  as  a Criterion  for  Solai  Proton  Event  Warning,"  AIAA  Joiiin.il. 
v 2.  p 2143.  1464 

6.  NASA  Program  Working  Papei . Preliminary  Warning  Critei  i.i  for  the  Solai  Panicle  A lei  t Netwoik,  b\ 

MD Lopez,  At  Bragg, and  JD Modisette,  196b 

7.  NASA  lechnical  Report  32-103.  Solar  Pioton  Forecast  System  and  Procedures  Used  Dining  the  Mai  met 
V Mission,  by  CC  Gonzalez  and  FI  Divita.  14(>8 

8.  Stiaka . RM.  and  Barron.  WR.  “Multifrei|uency  Solai  Radio  Bursts  as  Pi  edict  ors  for  Proton  I vents." 
AGARD  Conference  Proceedings  No  44.  l»(>4 

0 AFCRI  Space  Forecasting  Reseat  cli  Note  2.  1 lie  Use  ot  Solai  Radio  Buislsas  Piedictors  ol  Pioton 
Event  Magnitude,  by  RM  St i aka . |47() 

10.  ( room.  D1  . "Forecasting  the  Intensity  of  Solai  Proton  Events  from  the  I line  Chaiacteiistics  of  Solai 

Microwave  Bursts."  Solar  Phy  sics,  v |4.p  |7|.  |47| 

11  AFCRI  r echnic.il  Report  72-0543,  Fotecaslmg  Peak  Pioton  Flux  and  PC  A I vent  Magnitudes  Using 
‘Flash-Pliase'  Integrated  Radio  Bmst  I lux  Density . b\  Dl  Newell.  |4’2 

12  AFCRI  Technical  Report  73-00X6.  Solai  Radio  Activity  in  \ugiist  |472.b\  IP  Castelli.  WR  Ration, 
and  J Aarons.  |473 


TABU  1.  RFVIFW  OF  ML TIIODS  I OR  PR!  DlCTlNti  HU  1NTFNSITYOF 
SOLAR  PAR  I It'Ll  I VI  NTS  BASFI)  UPON  MIC  ROWAVI  BURST  DATA. 


Investigator 

Frequency 

Predictor  (Burst  Parameter) 

Predtctand 

Number 
of  1 veins 

Wohhor  ( 1 '>k.')'4 

10(111/ 

Fnergy 

Time-integrated 
intensity  of  >10- 
MeV  particles. 

15 

Fletcher 
( I'lM)5 

2.8  (111/ 

Duration,  peak  !lu\  density, 
rise  rate,  ratio  of  the  duration 
of  the  hurst  to  the  duration  of 
the  burst  at  one-half  peak 
flux  density. 

Tune-integrated 
flux  of  >30  MeV 
protons. 

I 

Lopez  et  al 
( 1 '•(!(')*' 

2.8  (ill/ 
.V75  till/ 

Fnergy.  energy  above  a 50-Sfu 
baseline,  energy  above  baselines 
of  10  and  20  percent  of  the 
burst's  peak  flux  density  . peak 
flux  density . duration,  mean 
flux  density . 

Time-integrated 
tin x of  >50  MeV 
protons. 

20 

(lon/ale/  and 
Divita  ( l%8) 

:.8(;ii/ 

Fnergy.  peak  flux  density, 
product  of  burst  duration  and 
peak  flux  density. 

Time-integrated 
tin  x of  >50-MeV 
protons. 

24 

Strata  and 
Barron 
( 1 '»0'))8 

500  Mil/  to 
8.5  (ill/ 

Total  integrated  radio-flux 
density  from  500  Mil/  to 
8.5  (ill/. 

Peak  proton  flux 
>10  MeV.  Peak 
riometer  absorp- 
tion at  50  Mil/. 

lb 

Stiaka  ( 1 ‘>70)° 

bOb  Mil/, 
1.4.  2.7. 

5.0,  S.S.  and 
15.4  (ill/ 

I nergy.  mean  flux  density  , 
peak  flux  density  . dotation. 

Peak  proton  flux 
>10  MeV  Peak 
riometet  absorp- 
tion at  50  Mil/. 

15 

('room  ( 1 ‘>7 1 )* 11 

5 to  20  (ill/ 

Mean  duration.* 

Peak  proton  flux 
>10  MeV. >50 
MeV.  '>>0  MeV 

57  1 

Newell!  1 ‘>7  2I1  1 

(>0(i  Mil/. 
1.4.  2.7. 

5.0,  8.8  (ill/ 

Semi-integiated  flux  density 

Peak  proton  flux 
>10  MeV.  Peak 
ttomelet  absorp- 
tion at  50  Mil/ 

lb 

t astelli  ol  al 

(l‘>75),: 

8.8  (ill/ 

Fnergy 

Peak  proton  Hit x 
''10  MeV 

21 

i 

('mom 4 11  considered  sevetal  single liei|tienc \ .mil  multilrei|iienc\  parameleis  ol  niictowave  butsts 
llo  concluded  ili.il  host  results  wcte  obtained  h\  using  tlto  effective  dotation  .is  a predictor 
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The  latter  of  these  is  an  extension  of  Strata’s  study  at  8.8  GHz  and  provides  the  primary 
input  to  a program  recently  developed  by  Smart  and  Shea  (in  preparation)  which  predicts 
the  time  behavior  of  the  Hare  particle  spectrum.  (A  third  method  of  predicting  the  intensity 
of  PCA  events  which  is  currently  in  use  at  AFGWC  was  developed  by  Kuek  (ref  1 3)  and  Kuck 
et  al  (ref  14).  It  is  based  on  the  parent-flare  soft  X-ray  emission.) 

The  main  difference  between  this  study  and  those  summarized  in  table  1 is  that  the 
data  sample  considered  here  is  much  larger.  First  of  all,  more  events  have  been  observed  since 
Castelli  et  al  (ref  I 2)  reported  their  results.  Second,  correlating  with  riometer  absorption 
rather  than  satellite-measured  particle  flux  allows  one  to  include  events  from  the  nineteenth 
solar  cycle.  As  an  added  bonus,  many  of  these  early  absorption  events  are  among  the  largest 
ever  observed.  The  predictor  used  is  that  which,  over  the  years,  has  emerged  as  the  most 
reliable  event  energy.  The  physical  reason  for  expecting  radio-event  energy  to  be  correlated 
with  peak  particle  flux  is  that  the  radio  energy  should  be  proportional  to  the  number  of  elec- 
trons accelerated  which  should,  in  turn,  according  to  neutrality  arguments,  be  approximately 
equal  to  the  number  of  positive  ions  which  are  accelerated.  Finally,  the  choice  of  the  radio 
frequency  for  this  study  was  dictated  both  by  Straka’s  preliminary  indication  that  low  micro- 
wave  frequencies  are  superior  for  predicting  peak  riometer  absorption  and  also  by  the  availa- 
bility of  the  excellent  2.8-GHz  (10.7-cm)  observations  by  Covington  at  Ottawa  throughout 
the  nineteenth  and  twentieth  solar  cycles. 


DATA  CONSIDERATIONS 

The  selection  criteria  for  the  events  were  as  follows: 

(a)  30-MHz  riometer  absorption  > 0.5  dB; 

(b)  Visible  hemisphere  parent  Hare;  and 

(c)  Parent  flare  observed  by  Ottawa  (nineteenth  solar-cycle  events  only). 

Criterion  (a)  was  imposed  because  0.5  dB  is  near  the  lower  limit  of  riometer  sensitivity. 
The  prediction  threshold  beyond  which  AFGWC  notifies  customers  of  an  impending  polar-cap 
event  is  1.0  dB  of  nighttime  absorption  and  2.0  dB  of  daytime  absorption. 

Criterion  (b)  is  necessary  because,  for  limb  events.  ( 1 ) radio  (and  optical)  source 
regions  are  occulted  by  the  disk  by  an  unknown  amount,  and  (2)  particle  propagation  effects 
may  play  a dominant  role. 

Criterion  (c)  was  adopted  because  of  the  availability  and  consistency  of  the  Ottawa 

data. 

By  the  mid  l%0s.  coverage  of  the  sun  at  or  near  2.8  GHz  was  virtually  complete  and 
all  twentieth  solar-cycle  events  which  met  criteria  (a)  ami  (b)  were  included  in  the  sample.  In 
all,  70  events  satislVd  the  criteria.  These  events  are  listed  in  tables  2 and  3.  For  all  events  in 
table  2,  the  radio  emission  of  the  parent  llare  met  or  exceeded  the  500-Sfu  peak  flux  density* 
“bell-ringer"  criterion  of  the  Air  Weather  Services  Space  Environmental  Support  System 
(SI  SS)  network  of  observatories  at  some  frequency,  f.  > 1415  MHz.  The  parent  flares  listed 
m table  3 did  not  produce  (or  would  not  have  produced)  a radio-event  warning. 


13.  A1WI  Uvlink.il  Note  \A  I R I II.  Prediction  of 'Polar  Cap  Absorption  (vents.  In  U A Ruck.  I'l’i 
14  AIWI  Uvltnica  I Report  .'l  I . Prediction  of  Polar  Cap  Ab  rrption  Events,  by  G A Kuck,  RR  Davis,  and 
(•J  Krause.  1 071 
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1 1 
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2B 
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TABLE  2.  (Continued). 


15  Apr  70 

23  Jul  70 

05  Nov  70 
1 1 Dec  70 

24  Jan  71 

06  Apr  7 1 
05  Mar  72 

28  May  72 

02  Aug  72 
04  Aug  72 

07  Aug  72 
30  Oct  72 

29  Apr  73 

03  Jul  74 

04  Jul  74 
10  Sep  74 
19  Sep  74 

05  Nov  74 
21  Aug  75 

30  Apr  76 


H-a  Optical 

Peak  Flare 

Time  Class 


2 .8-011/  30-MHz 

Me  Math  F.vent  Riometer  * 

Solar  Plage  Energy  Absorption  Observa- 

Coordinates  Region  (Sfumtin)  (dB)  lory  Notes 


N13.WS6 


N09.E09 

S12.E36 

N16.W02 

N18.W49 

S19.WS0 

S07.E43 

N09.E30 

N14.E28 

N 14. EOS 

N14.W37 

S10.W03 

N14.W73 

S14.E08 

S16.W08 

N10.E61 

N09.Wb2 

S12.W78 

N28.W76 

S09.W47 


14  1 10 


379  200 
1 24  630 


68  640 


34  695 
49  350 


*0  = Ottawa;  N = Nagoya;  M = Manila;  SG  = Sagamore  Hill;  Ne  = Nederhorst  . P = Penticton;  SL.  = Slough; 

L = La  Posta 

NOTES  TO  TABLE  2 

( 1 ) A riometer  absorption  value  of  7.8  dB  was  also  reported. 

(2)  The  midpoint  of  the  range  of  riometer  absorption  values  w'as  used  in  the  correlations. 

(3)  A sudden  commencement  of  a great  geomagnetic  storm  (GMS)  occurred  on  1 1 May  at  2330  IT.  near 
the  time  of  maximum  riometer  absorption. 

(4)  A sudden  commencement  on  17  July  at  1638  IT  which  preceded  a great  GMS  is  associated  with  the 
observed  particle  flux  maximum  of  this  flare.  Peak  riometer  absorption  occurred  on  1 7 July  at  ' 10  bouts  l 1 

(5)  A sudden  commencement  of  great  GMS  occurred  on  08  May  at  042 1 IT.  neat  the  peak  of  the  PGA 


16)  The  peak  time  listed  is  that  of  the  2.8-Gllz  emission. 

(7)  The  radio  observations  lot  ibis  event  ate  at  2.0  (ill/. 

(8)  A GSSC  occurred  on  29  Aug  al  1315  IT.  '-22  hours  at  let  the  Hate  and  - I 2 hours  helot  c l he  l’<  \ 
event  maximum. 


IABLI  2.  (Continual). 


('>)  Following  Svestk.i  and  Simon  (ref  15).  the  PCA  maximum  is  taken  at  2550  IT. 

( 10)  A GSSC  occurred  on  24  May  at  1 72b  L’  1,  " 22  hours  after  the  flare  and  " I')  hours  before  the  P('A 
event  maximum.  A second  SC  occurred  on  25  May  at  "12  hours  LIT.  near  the  peak  of  the  PC  A event. 

(11)  The  radio  observations  for  this  event  are  at  5.0  (111/. 

(12)  The  2 .S -Ci  11/  event  energy  was  measured  from  published  Manila  strip-chart  records. 

( 15)  The  Manila  strip-chart  record  was  not  available  for  this  event. 

( 14)  Geomagnetic  storm  modulation  with  sudden  commencement  on  05  Dec  at  Oo55  IT  is  listed  as  a deti- 
nue source  of  particles  for  a subsequent  PC  A event.  The  peak  riometer  absorption  for  the  Hare- 
associated  event  occurred  on  05  Dec  at  l)'>42  UT. 

( 1 5)  A major  magnetic  storm  began  at  05  hours  IT  and  had  its  maximum  at  IS  hours  IT.  near  the  peak  of 
the  PC A event . 

(lb)  A sudden  commencement  occurred  at  0*47  UT  on  lb  April,  near  the  peak  of  the  PC  A event. 

( I 7)  The  absorption  "spike"  of  4.7  dB  is  associated  with  a GSSC  at  2550  IT  on  24  July.  " 24  hours  after 
the  prompt  phase  of  the  particle  event.  An  earlier  sudden  commencement  which  occurred  on  24  July 
at  "II  hours  UT  was  probably  associated  with  a 21)  flare  ( N0l> .1- 52)  with  maximum  at  0052  IT  on 
22  July. 

(18)  A sudden  commencement  occurred  at  004b  IT  on  07  Nov,  about  b hours  before  the  peak  of  the  low- 
energy  particle  event . 

( I'D  The  radio  event  was  observed  shortly  before  sunset  at  Manila. 

(20)  Castelh  (GSDB.  Appendix  A)  and  Van  llollebeke  et  al  list  a -N  event  (S07.F40)  with  maximum  at 

1 141  IT  as  an  additional  source  of  particles  for  this  event.  A sudden  commencement  was  avoided  at 
"21  hours  IT  on  0b  March.  I hour  before  the  PC  A maximum  was  recorded  at  Thule  and  "57  hours 
after  the  flash  phase  of  the  IB  flare. 

(21)  A sudden  commencement  associated  with  an  earlier  flare  on  02  Aug  ( 1 B;  N 14.1-55 ; 04 10  II  I ) occurred 
at  01  l'>  UT  on  04  Aug.  The  sudden  commencement  associated  with  this  flare  occurred  on  04  Aug  at 
0221  IT  (ref  lb).  The  50-Mll/  absorption  reached  a peak  of  8.5  dB  al  ' 0b  hours  UT  on  04  Aug.  Fol- 
lowing Castelh  and  Van  llollebeke  et  al,  the  1 B flare  on  02  Aug  is  considered  to  be  a contributing 
source  of  particles  for  the  PC  A event. 

(22)  The  Thule  50-Mll/  riometei  saturated  at  "12  hours  IT. 

(25)  Sudden  commencements  occurred  on  08  Aug  at  2554  UT  and  on  ()'>  Aug  at  0057  UT.  l ire  first  SC  is 
associated  with  this  flare.  Maximum  riometei  absorption  was  recorded  at  Thule  at  " 22  hours  on 
08  Aug. 

(24)  Fhe  flare  association  for  this  event  is  after  Castelh.  Van  llollebeke  et  al  identify  an  earlier  Date  (2N, 

SI 0.1- 05 ) repor ted  by  lluancayo  as  the  particle  source.  This  flare  began  al  1 544  UT  on  2'*  Oct  and 
lasted  over  4.5  hours.  Radio  activity . during  this  time  period,  consisted  of  a gradual  rise  and  fall  event 
beginning  at  I 55(1  UT  and  ending  at  2 1 08  UT  with  superimposed  complex  events  peaking  at  I b!  5 l T 
and  1755  IT.  In  all.  the  2.8-011/ energy  totaled  nearly  8000  Sfu-min.  Caslelli  lists  this  eatliei  flare  as 
a possible  contributing  source  of  the  panicle  event.  The  flare  in  the  table  was  chosen  as  the  principal 
particle  source  because  il  is  a mote  likelv  candidate  to  have  produced  the  shock  icsponsible  lot  the  SSC 
on  51  Octal  Ib54  UT.neat  the  peak  of  the  absoiption  event.  The  tadto  emission  of  this  flare  was 
moderately  intense  at  the  liighei  microwave  frequencies  (Sp  peak  flux  density  '>40  Stu  at 
f - 8.8  (ill/).  Also  lluancay  o was  the  onlv  observatory  of  t lu ee  on  patrol  between  I 5 horns  l 1 and 
20  hours  IT  on  2'*  Oct  which  assigned  an  importance  classification  of  2N  to  the  cent i al  disk  activity 
I'he  oilier  two  stations  reported  only  subflares. 

( 25  ) Fins  flare  was  the  source  ol  etiei get ic  particles  w Inch  produced  a ground -level  event  but  a p patent  l\  it 
produced  relatively  lew  low-euetgy  protons. 
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FAB  LI'  2.  (Continued). 

(7b)  Following  Caste  III.  possible  contributing  Hares  had  maxima  at  071  b CT  (IN.  S13.E  17)  on  07  July  and 
0840  CT  (71).  S14.E08)  on  07'  July. 

(77)  Castelli  lists  a total  of  seven  events  occurring  on  04  and  05  July  which  may  have  contributed  to  this 
particle  event.  Of  these,  the  event  in  the  table  was  the  most  energetic  at  centimetre  wavelengths.  Sud- 
den commencements  occurred  on  05  July  at  1450  UT  and  on  0b  July  at  0577  CT.  It  is  believed  that 
the  first  SC  is  associated  with  the  previous  flare  in  the  table  and  the  second  with  the  Hare  under  discus- 
sion. The  PCA  maximum  was  recorded  on  05  July  at  70  hours  CT.  It  is  suggested  that  the  magnetic 
field  of  the  shock  front  responsible  for  the  first  sudden  commencement  inhibited  the  prompt  low- 
energy'  particle  response  of  this  flare. 

15.  Svestka,  Z,and  Simon,  P,  "Proton  Flare  Project."  Solar  Physics,  v 10,  p 3,  14b4 

I b.  Dryer,  M,  E via  tar.  A,  Frohlich,  A.  Jacobs.  A,  Joseph,  JH,  and  Weber.  EJ . “Interplanetary  Shock  Waves 
and  Comet  Brightness  Fluctuations  During  June-August  1477,"  Journal  of  Geophysical  Research, 
v 80,  p 700 1 . 1 475 


TABLE:  3.  EVENTS  FOR  WHICH  THE  PARENT-FLARE  PEAK  RADIO  EMISSION 

<500  Sfu(f  >1415  MHz). 


TABLK  3.  (Continued). 


NOTES  TO  TABLt:  3 

( 1 ) The  nometer  absorption  was  superimposed  on  an  earlier  absorption  event. 

(2)  Geomagnetic  storm  modulation,  with  sudden  commencement  at  2348  UT  on  1 5 Feb.  is  listed  as  a 
definite  contributing  source  of  particles. 

(3)  The  radio  observations  are  at  2.0  (ill/. 

(4)  The  Manila  strip-chart  record  was  not  available  for  this  event. 

(5)  The  peak  time  is  that  of  the  2.8-GII/  emission. 

(6)  Observations  at  Manila  were  terminated  by  sunset. 

(7)  Radio  observations  are  at  3.0  GH/.  The  peak  time  is  that  of  the  3.0-GHz  emission. 

(8)  Van  Hollebeke  et  al  list  a I B Hare  (SI4.K48)  with  maximum  at  1128  UT  as  the  source  of  particles  for 
this  event.  However,  the  flare  in  the  table  is  a more  likely  candidate  to  have  produced  the  shock  front 
responsible  for  the  sudden  commencement  on  08  Mar  at  ~15  hours  UT  that  was  associated  with  the 
peak  of  the  riometer  absorption  event.  Neither  event  was  very  impressive  at  radio  wavelengths. 

(0)  The  radio  observations  are  from  Toyokawa  at  2.0  GHz. 

(10)  A sudden  commencement  was  observed  on  2 1 Jan  at  1151  UT,  near  the  peak  of  the  PCA  event.  0 

(11)  The  flare  association  for  this  event  is  after  Castelli. 

(12)  Coordinated  observations  of  these  flares  during  the  C'lNOF  program  ( ref  1 7-20)  suggest  that  both 
flares  may  have  contributed  to  the  particle  event.  In  Figure  1 . the  sum  of  the  event  energies  is  used. 

Sudden  commencements  occurred  at  0630  UT  and  1311  UT  on  I 7 June.  The  second  SC  has  been 
associated  with  the  flare  at  0958  UT  on  15  June.  The  Shepherd  Bay  riometer  peaked  at  18  hours  UT 
on  1 7 June. 

I 7.  Dodson.  HW.and  Hedeman,  RF.  “A  Small  Subflare  in  the  CINOF  Program:  June  1972.  I4l*.  |9*\  3(>m 
UT.”  AFCRI.  Technical  Report  75-0437.  1975 

18.  Bruzek.  A,  “The  June  1972  Event."  AFCRI  Technical  Report  75-0437.  1975 

1 9.  Shea.  MA.  and  Smart . DF.  “An  Analysis  of  Solar  Particle  Data  Obtained  During  CINOI  AFCR I 
Technical  Report  75-0437.  Results  Obtained  During  the  Campaign  for  Integrated  Observations  of 
Solar  Flares,  1975 

20.  Mihalov.  JD.  and  Wolfe.  .III.  “Pioneer  10,  (■>,  and  9 Solar  Plasma  Observations  During  5-29  June  1972." 

AFCRI,  Technical  Report  74-027 1 , Compilation  of  Solar  Particle  and  Interplanetary  Measurements 
Acquired  During  I he  Campaign  lot  Integrated  Observations  of  Solai  Flares.  1974 

The  2.8-Gllz  energy  listed  in  tables  2 and  3 is  that  of  the  main  radio  event;  ie.  it  does 
not  include  the  energy  of  the  precursor  and  post  burst  increase  which  often  accompany  large 
bursts.  The  energies  have  been  calculated  from  the  burst  durations  and  mean  flux  densities 
listed  in  Solar  Geophysical  Data  ( S(  il)l.  the  Quarterly  Bulletin  of  Solar  Activity  (QBSA ).  and 
the  Geophysics  and  Space  Data  Bulletin  (GSDB).  As  pointed  out  by  Straka  and  Barron  (ref  8), 
using  the  relationship 

I = SM  • T ( 1 1 

where 

I Burst  Energy  (Sfirmin). 

Sm  = Mean  Flux  Density  (Slu), 

I = Duration  ( min ). 


It) 


I 


to  determine  event  energy  is  not  without  its  risks.  The  eyeball  methods  often  used  to  obtain 
the  mean  flux  densities  of  complex  bursts  can  result  in  significant  errors  as  much  as  a factor 
of  two.  For  36  of  the  70  events,  the  2.8-Gllz  energy  was  calculated  from  the  accurate  Ottawa 
and  Penticton  burst  parameters.  In  addition,  profiles  of  large  bursts  are  often  published  in 
data  compilations,  enabling  us  to  measure  directly  the  area  under  the  curve  for  these  events 
and  further  reduce  probable  error.  The  energies  of  I 2 more  events  were  determined  in  this 
manner.  Following  Wefer  (ref  21  ),  the  energy  is  given  in  the  convenient  units  of  Sfu'tnin, 
rather  than  in  the  more  conventional  units  of  joules  per  square  metre  per  hertz  (J  • • II/-' ). 

As  a final  comment  on  the  radio  data,  it  should  be  noted  that  Ottawa  and  La  Posta  observe 
the  sun  at  2.8  (ill/,  while  the  nominal  frequency  of  observation  at  Sagamore  Hill  and  Manila 
is  2695  MHz.  In  actual  practice,  of  course,  this  slight  difference  in  frequency  is  generally  not 
significant  and  in  this  report  all  such  observations  are  referred  to  as  being  at  2.8  GHz.  For  a 
few  twentieth  solar-cycle  events,  observations  were  not  available  at  2.8  GHz;  for  these  cases, 
the  parameters  of  the  next  closest  reported  frequency  were  used  to  calculate  the  energy.  All 
such  exceptions  are  noted. 

For  events  occurring  before  1 970,  the  peak  riometer  absorption  value  and  the  identi- 
fication of  the  parent  Hare  were  taken  from  reference  22.  Only  those  catalog  events  with 
definite  parent-flare  associations  were  included  in  the  sample.  For  PCA  events  occurring 
since  1969  for  which  there  was  some  question  about  the  assignment  of  a parent  flare,  the  list 
of  Van  Hollebeke  et  al  (ref  23)  and  appendix  A of  the  GSDB  were  consulted.  For  these  events 
the  source  of  the  parent-flare  identification  is  indicated  in  the  notes.  Riometer  absorption 
values  for  the  post-1969  events  were  taken  primarily  from  the  GSDB  and  also  from  Cormier’s 
report  on  events  observed  at  Thule  (ref  24). 

Large,  favorably  located  (45° E <2  vV  ^ 90°W)  flares  often  produce  a delayed  increase 
in  riometer  absorption  associated  with  a geomagnetic  storm  sudden  commencement.  The 
sudden  commencement  is  caused  by  the  flare-produced  shock  front  hitting  the  earth’s  mag- 
netic field.  Low-energy  protons,  trapped  in  the  magnetic  field  of  the  shock,  account  for  the 
enhanced  absorption.  Typical  delay  times  between  flare  maximum  and  geomagnetic  storm 
onset  range  from  24  to  48  hours  for  an  event  at  45°W  and  from  48  to  72  hours  for  a flare  at 
45°E.  There  is  considerable  scatter  about  these  time  ranges  depending  upon  both  the  size  of 
the  flare  and  the  eondition  of  the  preflare  interplanetary  medium.  Thus,  identifying  the  flare 
responsible  for  a delayed  absorption  increase  is  inherently  more  difficult  than  identifying  the 
source  of  a prompt  particle  event.  This  is  an  important  consideration  since  the  delayed  riome- 
ter absorption  peak  associated  with  the  sudden  commencement  often  exceeds  the  earlier 
prompt-particle  induced  enhancement.  The  times  of  sudden  commencements  which  may 
have  influenced  or  determined  the  time  of  the  PCA  maximum  are  given  in  the  notes  to  tables 
2 and  3.  These  times  are  taken  from  reference  22  for  events  occurring  before  1 970  and  from 
the  SGI)  Prompt  Reports  or  reference  25  for  later  events.  Unless  it  is  otherwise  noted  or  is 
obvious  from  the  delay  times  that  it  cannot  be  so.  il  is  lo  be  assumed  that  the  sudden  com- 
mencement is  associated  with  the  flare  in  the  table. 

21 . PSD  Department  of  Astronomy  Scientific  Report  02<>.  A Statistical  Study  of  5480  Microwave  Sol.ii 
Radio  Events,  bv  M Wefer,  1973 

22.  Svestka.  Z. and  Simon,  P.  ed,  “Catalog  of  Solar  Paiticle  Events.  I '155-1 0(»9,"  Astrophysics  and  Space 
Science  I ibrary . I)  Keidel  Publishing  Company . Dordrecht.  Holland.  Il>75 

23.  Van  lloleheke,  MAI.  Ma  Sung.  I S.  and  McDonald.  I It,  "The  Variation  of  Sol.u  Proton  Energy  Spectra 
and  Si/e  Distribution  with  lleliolongitude. " Solai  Physics,  \ 41,  p |S‘>,  I ‘>75 

24.  AE’CKI  l echnic.il  Report  73-(K)b().  Elude  Kionietei  Observations  ol  Polai  Cap  Absorption  Events 
( 10(0-1972).  In  R.I  Cornuer . 1 11  '3 

25.  Van  I lollebeke.  MAI.  Wang.  I K . and  McDonald,  I H.  "\  Catalog  ol  Solai  t osmu  Ra\  I veins.  IMPS  l\ 
and  V (May  PR, 7-Dec  I *>72 " < • I SC  \ (.(vL74  27.  1074 
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ANALYSIS.  RESULTS,  DISCUSSION 


As  a first  approach  to  the  problem,  peak  riometer  absorption  (A)  was  plotted  against 
2. 8-till/  event  energy  (L)  for  the  70  events  listed  in  tables  2 and  2.  A power-law  relationship 
between  the  variables  was  obtained  by  the  standard  least -squares  technique.  I he  equation  ol 
the  line  in  figure  I is 

A (dB)  = 0.0040  EiSfu-min)0-3498 

(2h 

R = 0.62 

According  to  this  prescription,  a 2. 8-C.lt/  burst  with  L ' 6000  Sfirmin  will  be  associated 
with  a principal  riometer  absorption  event.  However,  inspection  of  a list  of  radio  events  com- 
piled in  SC’.D  or  elsewhere  reveals  many  such  events  which  do  not  have  proton  association. 

The  500  Sfu  SUSS  "bell-ringer”  criterion  is  designed  to  lower  the  PC  A alert  false-alarm  rate 
by  calling  attention  to  only  those  Hares  which  can,  on  the  basis  of  statistics,  be  expected  to 
produce  protons,  f igure  2 is  a plot  of  A versus  E for  the  55  events  in  table  2.  The  best  power- 
law  fit  to  the  data  is  given  by 


A (dB)  = 0.0709  E(Sfu-min)0J70: 
R =0.61 


(5) 


This  equation  is  not  significantly  different  from  equation  (2).  In  order  to  achieve  an  accept- 
able false-alarm  rate  it  is  necessary  to  take  the  process  a step  further.  In  figure  2.  22  of  the 
24  (94  percent)  principal  absorption  events  are  related  to  Hares  which  had  2. 8-till/  energy  in 
excess  of  2000  Sfu  * min.  If  we  restrict  the  sample  to  include  only  those  events  with  parent 
Hare  2.8-C.ll/  energy  > 2000  Sfu*  min.  we  have  figure  2.  The  equation  of  the  solid  line  is 
given  by 


A (dB)  = 0.01  16  E(  Sfu  -min)0*5555 
R = 0.59 


(4) 


I bis  formula  provides  greater  discrimination  between  those  Hares  which  are  capable  of  caus- 
ing a principal  polar-cap  event  to  occur  and  those  which  are  not.  It  predicts  that  protons 
ejected  from  Hares  with  E_>  I I 000  Sfu  • min  at  2.8  (ill/  will  produce  > 2.0  dB  of  20-Mil/ 
riometer  absorption.  As  a comparison,  in  appendix  A of  the  (iSI)B.  C'astelli  adopts  a burst 
energy  of  It)-'  .1  * in--  *11/“'.  or  I 667  Sfu  * min  at  a given  microwave  frequency  as  one  of 

the  criteria  which  a flare  should  meet  to  be  considered  as  a possible  source  of  protons  (at  least 
1 proton  * cm--  • s-'  • sir-*  ) for  E > 1 0 MeV.  or  after  Kuck  el  al  ( ref  1 4 1 A % 0.6  d B.  In 
table  4.  the  2. 8-dll/  energy  threshold  and  the  modified  U-shaped  spectral  criterion  (ref  2) 
are  compared  as  yes  or  no  predictors  of  principal  IVA  events  for  the  1968-1969  time  frame. 

I Ins  table  contains  a complete  list  for  these  2 years  of  relatively  high  solar  activity  of  all 
2. 8-till/  bursts  with  E * I 0 1)00  Sfu  • min.  all  radio  events  which  satisfied  the  modified  U cri- 
terion. and  all  > 0.5-dB  riometer  absorption  events  which  had  definite  parent-flare  associa- 
tions.* for  the  period  in  question,  it  can  be  seen  that  the  2. 8-till/  energy  threshold  compares 
reasonably  well  as  a yes-no  indicator. 


I lie  false  .il.iitn  i.nlio  events  in  i lie  table  were  taken  limn  observations  ol  i lie  Sagamore  Hill  ami  Manila 
observatories  published  in  (lie  ( .SDH  Ir  is  possible  tliai  a lew  sm li  events  inclined  dining  gaps  in  the  eov 
eiage  ol  these  two  nhseiv atones  and  thus  weie  missed  \s  m tables  and  ' whenever  possible . the  Ottawa 
bin  si  p.n  a met  et  s weie  used  to  c.  denial  e the  eneigy  ot  il  was  me.isliied  diiectlv  I tom  I lie  sti  ip-vlini  t ivcoi  ds 


li 


2.8  GHz  EVENT  ENERGY  (Sfu-min) 

Figure  3.  Riometer  absorption  versus  event  energy  (table  2 events 
with  F.  > 3000  Sfu  ■ min). 

The  dashed  lines  in  figure  3 indicate  events  which  lie  within  factors  of  2 and  3 of  the 
least-squares  line.  Since  the  10-MeV  proton  flux  varies  as  the  square  of  the  riometer  absorp- 
tion (ret  2b),  these  factors  correspond  to  factors  of  4 and  0 in  the  low-energy  proton  flux. 

Of  the  32  principal  PC  A events  in  the  figure,  22  (00  percent)  would  have  been  predicted 
within  a factor  of  2 anil  28  (88  percent)  would  have  been  called  within  a factor  of  3. 

In  conclusion,  it  has  been  shown  that  the  use  of  equation  (4).  in  conjunction  with  the 
modified  U-shaped  spectral  (C'astelli)  criterion,  can  be  expected  to  provide  an  accurate  esti- 
mate of  the  si/e  of  an  impending  principal  absorption  event  with  a minimum  of  false  alarms. 
A few  comments  on  the  limitations  of  this  study  are  in  order.  First  of  all.  no  propagation 
correction  has  been  built  into  this  prediction  tool.  In  figures  1-3.  the  eastern  hemisphere 
events  are  scattered  uniformly  about  the  respective  least -squares  lines;  however,  there  is  a 
marked  east-west  asymmetry  in  the  distribution  of  the  70  events  on  the  solar  disk  ( fig  4) 
Second,  before  any  forecasting  method  can  be  accepted  as  valid,  it  must  pass  an  operational 
test.  The  Polar  Cap  Absorption  events  of  the  twenty-first  solar  cycle  will  provide  the  neces- 
sary independent  data  sample. 


3(i.  Juil.iv . Kl).  ami  Adams.  ( >\V  . "Riometei  Measuieinenis.  SoI.ii  Pinion  Intensities  and  Radiation  Dose 
Rates."  Planet . Space  Science.  \ 17,  p 1313.  I'bet 


TABLE  4.  PREDICTION  METHOD  APPLIED  TO  1968-1969  EVENTS. 


EAST  WEST 

SOLAR  LONGITUDE 


Figure  4.  Longitudinal  distribution  of  parent  Hares. 
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APPENDIX  A:  PARTICLE  PROPAGATION  CORRECTION 

The  proton  prediction  program  developed  by  Smart  and  Shea  for  AFGWC  contains  a 
correction  tor  the  effect  of  the  interplanetary  magnetic  field  on  particle  propagation,  finis 
equation  (4)  must  be  modified  if  it  is  to  be  used  as  an  input  to  their  program;  otherwise  the 
predictions  made  will  generally  be  too  low. 

A given  proton  Hare  will  produce  a maximum  particle  flux  (measured  at  earth)  when 
it  is  located  at  the  solar  foot  point  of  the  Archimedes  spiral  which  intersects  the  earth,  and 
less  when  located  elsewhere.  According  to  Smart  and  Shea’s  prescription,  the  reduction  in 
flux  with  distance  of  the  flare  from  this  footpoint  goes  as  e-^,  where  t)  is  the  angle  in  radians 
between  the  heliographic  longitude  of  the  Hare  region  and  the  longitude  of  the  footpoint. 
both  taken  at  the  time  of  peak  riometer  absorption.  The  solar  longitude  of  the  footpoint  of 
the  Archimedes  spiral  field  line  which  connects  to  the  earth  is  given  by 

404.1 

Solar  longitude  ot  tootpoint  = — — — (radians)  ( A 1 ) 

vs 

where  Vs  is  the  solar  wind  speed  in  ktirs-' . The  National  Space  Science  Data  Center  has 
compiled  a comprehensive  list  of  the  solar  wind  measurements  available  for  the  period  from 
1964  to  1973.  For  events  occurring  before  1964  (or  since  during  data  gaps),  a nominal  Vs  of 
404.1  km  • s-1  may  be  used,  placing  the  footpoint  at  1 radian  or  ' 57°W.  Smart  and  Shea’s 
program  takes  into  account  other,  higher  order  corrections,  such  as  the  effect  of  Hare  latitude, 
but  the  longitude  effect  is  clearly  dominant. 

The  procedure  for  incorporating  the  propagation  correction  into  the  basic  result  of 
this  paper  (or  any  other  forecast  scheme)  so  that  it  can  be  utilized  in  AFGWC’s  proton  predic- 
tion program  is  as  follows; 

A.  For  each  event  in  the  data  base: 

1 . Convert  riometer  absorption  (A)  into  >5.2-MeV  proton  flux  (J)  using 


J (protons  • cm 


-1 


sr 


-1 


A(dB)2 
7T(0. 1 1 5)i 


(a:i 


after  Stroseio  and  Sellers  (ref  27).  (The  >5.2-MeV  peak  proton  flux  is  used 
in  these  correlations  because  of  its  relative  independence  of  the  slope  of  the 
proton  energy  spectrum  (ref  27.) 


Determine  0 at  the  time  of  peak  riometer  absorption,  taking  solar  rotation 
into  account. 


Multiply  the  >5.2-MeV  flux  measured  at  the  earth  by  e-',)  to  obtain  the 
flux  at  1 au  out  on  the  Archimedes  spiral  which  connects  to  the  flare  region 


B.  Determine  the  relationship  between  the  2.S-GI1/  event  energy  (1  )(or  any  other 
parameter)  and  the  particle  flux  from  (A)  by  the  least-squares  method. 

The  result  of  ( B)  for  the  47  events  considered  in  this  paper  is  given  in  figure  \ I and 
equation  ( A3). 

11.  AFCRl.T  echnival  Report  75-()4<e>,  Hie  Calculation  of  Rioiuctci  Absoiption  ami  an  Xppiosinutc  l on 
neclion  Between  Kiometei  Absorption  ami  Solat  Proton  I luxes  Dining  Night-time  IV  A I vents.  In  St  \ 
Stroseio  ami  It  Sellers.  I '>75 
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rile  improvement  in  the  correlation  which  one  would  hope  to  see.  as  affirmation  of  the  valid- 
it>  of  the  propagation  correction,  is  not  apparent.  It  should  he  noted  that  Newell  (ref  I I ) 
found  a similar  lack  of  improvement  when  he  applied  the  propagation  correction  of  Kuck  et  al 
(ret  14)  to  his  list  of  events.  Smart  and  Shea  kindly  provided  a test  run  of  their  program  for 
the  47  events  using  equation  (A3)  as  the  initial  input.  The  predictions  and  observations  are 
compared  in  figure  A2.  Twenty-nine  events,  or  (■> 2 percent,  lie  within  a factor  of  2 of  the 
optimum  line.  Those  events  falling  outside  the  acceptable  range  tend  to  be  overpredicted.  It 
is  perhaps  significant  that  only  one  of  the  I 7 events  with  eastern  hemisphere  parent  flares,  ie, 
those  events  to  which  the  largest  correction  factors  are  applied,  is  overpredicted  In'  this  method. 


Figure  Al . >5.2-MeV  peak  proton  flux  at  I au  on  the  Archimedes  Figure  A2.  Results  of  the  Smart  and  Shea  prediction  program 

Spiral  with  footpomt  in  the  parent-flare  region  versus  event  energy.  using  equation  (A3)  as  the  initial  input. 
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